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SUMMARY 


The  method  developed  in  NACA  TN  No.  995  has  "been  slightly  modified 
and  extended  to  include  flows  with  circulation.  The  essential  feature  of 
the  modified  method  is  that  in  analytic  continuation  of  the  solution  the 
alteration  of  the  singularities  of  the  incompressible  solution  due  to  the 
presence  of  the  hypergeometric  functions  has  been  taken  into  account. 

It  was  found  that  for  finite  Mach  number  the  only  case  in  which  the 
nature  of  the  singularity  of  the  incompressible  solution  can  remain 
unchanged  is  for  a  ratio  of  specific  heats  equal  to  -1. 

Two  particular  flows,  one  having  a  finite  circulation  and  the  other 
having  zero  circulation,  have  been  studied.  Both  flows  were  derived  from 
the  incompressible  flow  about  an  elliptic  cylinder  of  thickness  ratio  0.60. 
The  free— stream  Mach  number  for  both  cases  was  taken  to  be  0.60  in 
order  to  avoid  the  appearance  of  limiting  lines.  The  pressure  distribution 
for  the  flow  without  circulation  has  been  compared  with  that  of 
incompressible  flow  over  approximately  the  same  body.  The  discrepancies 
between  the  exact  results  and  those  predicted  by  the  approximate  Von  Karman- 
Tsien  and  G-lauert— Prandtl  formulas  are  so  wide  as  to  show  definitely  that 
in  this  case  the  effect  of  geometry  cannot  be  ignored,  as  is  done  in  both 
approximate  formulas.  In  general,  it  seems  that  the  effect  of  geometry 
cannot  be  neglected  and  the  conventional  "pressure-correction"  formulas 
are  not  valid,  even  in  the  subsonic  region  if  the  body  is  thick,  especially 
if  there  is  a  supersonic  region  in  the  flow. 


INTRODUCTION 


This  report  is  a  continuation  of  the  work  reported  in  NACA  TN 
No.  995.  The  method  developed  in  that  report  has  been  slightly  modified 
and  extended  in  the  present  report  to  include  flows  with  circulation. 

The  general  concept  and  method  is  outlined  without  the  mathematical 
details  in  part  I.  The  essential  feature  of  the  modified  method  is 
that  in  analytic  continuation  of  the  solution  the  alteration  of  the 
singularities  of  the  incompressible  solution  due  to  the  presence  of  the 
hypergeometric  functions  has  been  taken  into  account,  as  fully  discussed 
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in  part  II.  It  was  found  that  foi  finite  Mach  number  the  only  case  in 
which  the  nature  of  the  singularity  of  the  incompressible  solution  can 
remain  unchanged  is  for  a  ratio  of  specific  heats  equal  to  —1.  Part  III, 
which  contains  a  discussion  of  the  improvement  of  convergence  of  the 
power  series,  contains  no  essentially  new  material  and  is  included  primarily 
for  the  sake  of  completeness.  Detailed  proofs  are  given  in  appendixes  A 
to  F. 


Part  IV  contains  the  study  of  two  particular  flows,  one  having  a 
finite  circulation  and  the  other  having  zero  circulation.  Both  are  derived 
from  the  incompressible  flow  about  an  elliptic  cylinder  of  thickness 
ratio  0.60.  The  free-stream  Mach  number  for  both  cases  is  taken  to  be 
0.60  in  order  to  avoid  the  appearance  of  limiting  lines.  The  pressure 
distribution  for  the  flow  without  circulation  has  been  compared  with  that 
of  incompressible  flow  over  approximately  the  same  body.  The  discrepancies 
between  the  exact  results  and  those  predicted  by  the  approximate 
Von  Karmdn—Tsien  and  Glauert-Prandtl  formulas  are  so  wide  as  to  show 
definitely  that  In  this  case  the  effect  of  geometry  cannot  be  ignored,  as 
Is  done  in  both  approximate  formulas.  In  general,  it  seems  that  the 
effect  of  geometry  cannot  be  neglected,  and  the  conventional  "pressure- 
correction”  formulas  are  not  valid,  even  in  the  subsonic  region  if  the 
body  is  thick,  especially  If  there  is  a  supersonic  region  in  the  flow. 

The  importance  of  this  result  cannot  be  overemphasized,  as  there  is  a 
widespread  tendency  in  engineering  practice  to  use  simple  pressure- 
correction  formulas  indiscriminately. 

This  work  was  conducted  at  the  Guggenheim  Aeronautical  Laboratory  of 
the  California  Institite  of  Technology  under  the  sponsorship  and  with  the 
financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics. 


I  -  CONCEPTS  AKD  METHODS 
General  Consideration  of  Transonic  Flows 


The  flow  of  a  compressible  Ideal  fluid  about  an  infinite  cylindrical 
body,  unlike  that  of  an  incompressible  fluid,  depends  on,  among  other 
conditions,  the  speed  or  Mach  number  at  infinity.  If  the  free— stream 
Mach  number  is  below  a  certain  value,  the  flow  pattern  will  be  very 
similar  to  that  of  an  incompressible  fluid  even  though  part  of  the  flow 
may  be  supersonic.  However,  as  soon  as  the  limiting  Mach  number  is 
reached,  the  situation  is  entirely  different.  The  phenomenon  of  major 
importance  in  this  new  situation  is  the  appearance  in  the  calculations 
of  limiting  lines  in  the  supersonic  region,  characterized  by  the  fact 
that  the  fluid  particles  there  experience  an  infinite  pressure  gradient. 

It  can  be  shown  that  if  the  assumptions  of  isotropy  and  of  irrotationality 
of  the  flow  are  not  rejected,  it  is  impossible  to  continue  the  solution 
beyond  these  singular  lines  (reference  1).  The  failure  of  potential  flow 
can  be  attributed  to  the  effects  of  viscosity  and  conductivity  of  the 
fluid.  Although  the  exact  relation  between  the  limiting  line  and  shock 
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wave  is  still  not  established,  there  is  reason  to  believe  that  in  most 
cases  the  theoretical  appearance  of  a  limiting  line  necessarily  implies 
the  existence  of  a  shock  wave.  Therefore,  in  practice  the  knowledge  of 
the  conditions  under  which  the  limiting  lines  appear  for  a  given  body 
is  of  paramount  importance. 


Chaplygin  Holograph  Method 

To  solve  the  problem  of  a  potential  flow  of  compressible  fluid, 
the  use  of  the  hodograph  method  was  first  suggested  by  P.  Moleribroek 
(reference  2)  and  later  by  S.  A.  Chaplygin  (reference  3).  The  advantage 
of  this  method  is  that  in  the  case  of  two-dimensional  potential  flow  it 
leads  to  a  linear  partial  differential  equation  instead  of  a  quasi-linear 
one,  such  as  obtained  in  the  physical  plane.  The  particular  solutions 
of  this  partial  differential  equation  were  found  to  be  very  similar  to 
those  for  an  incompressible  fluid  except  that  each  consists  of  a  hyper- 
geometric  function  containing  the  Mach  number  as  a  parameter.  Since  the 
equation  is  linear,  a  general  solution  can  be  constructed  for  any  assigned 
domain  by  the  principle  of  superposition.  The  difficulty  in  connection 
with  this  method  is  that  the  solution  so  obtained  in  the  hodograph  plane 
may  not  transform  to  a  "good”  aerodynamic  body.  The  problem  of  con¬ 
structing  a  solution  in  the  hodograph  plane  which  corresponds  to  a 
desired  body  in  the  physical  plane  has  proved  to  be  extremely  difficult. 

The  difficulty  was  partly  solved  by  constructing  a  solution  which, 
in  the  limit  of  zero  Mach  number,  reduces  to  a  known  incompressible  flow. 
Using  this  idea,  Chaplygin  studied  the  subsonic  motion  of  a  gas  jet  and 
F.  Ringleb  (reference  4)  calculated  the  flow  around  a  sharp  edge. 

Following  the  same  principle,  Tsien  (reference  5)  for  the  first  time  solved 
the  problem  of  a  subsonic  flow  about  a  closed  body,  by  simplifying  the 
differential  equation  in  such  a  manner  that  the  difference  between  the 
compressible  and  incompressible  solutions  appears  only  in  a  modification 
of  the  speed  scale.  This  method  was  later  generalized  by  L.  Bers 
(references  6  and  7)  to  include  the  flow  with  circulation.  In  the  general 
case  where  the  ratio  of  specific  heats  y  is  not  equal  to  — 1,  this 
simple  form,  of  speed  distortion  does  not  exist;  therefore  the  method 
must  be  extended. 


Critique  of  a  Previous  Method  by  Tsien  and  Kuo 

In  reference  8  the  flow  past  a  closed  body  as  well  as  over  a  wavy 
surface  of  a  gas  with  a  characteristic  constant  y  greater  than  unity 
was  considered.  The  method  was  based  on  the  determination  of  two 
functions:  the  stream  function  and  the  transformed  potential  introduced 
through  a  Legendre  transformation.  When  these  two  functions  are  properly 
chosen,  the  coordinate  functions  x(q,8)  and  y(q,0)  are  given  by 
differentiation.  Consequently,  with  the  aid  of  the  stream  function,  the 
flow  pattern  can  be  calculated. 
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The  difficulty  of  the  scheme  is  twofold.  First,  whenever  a  stream 
function  and  a  transformed  potential  satisfying  the  conditions  of  conti¬ 
nuity  are  constructed,  they  are  required  to  satisfy  the  further  condition 
of  "compatibility"  to  ensure  that  they  actually  do  belong  to  the  same 
flow.  In  the  case  of  circulation— free  flow,  this  condition  brings  no 
difficulty.  However,  for  a  circulatory  flow  there  seems  to  be  an  insuf¬ 
ficient  number  of  arbitrary  constants  to  meet  this  stringent  requirement. 
Second,  the  coexistence  of  two  indirectly  connected  functions,  such  as 
the  stream  function  and  the  transformed  potential,  automatically  intro¬ 
duces  two  groups  of  hypergeometric  functions  and  doubles  the  number  of 
sets  of  unknown  coefficients  defining  the  power  series.  All  these  compli¬ 
cations  result  in  a  large  amount  of  labor. 

In  order  to  simplify  the  procedure,  the  velocity  potential  will  be 
Introduced  instead  of  the  transformed  potential.  Inasmuch  as  both  the 
stream  function  and  the  velocity  potential  are  directly  connected  by 
simple  differential  equations,  the  determination  of  one  leads  uniquely  to 
that  of  the  other.  The  difference  between  the  particular  integrals  of 
these  functions  is  simply  as  follows.  For  the  stream  function  the  par¬ 
ticular  integral  contains,  aside  from  the  trigonometric  function,  the 

factor  qvF^  av,bv  jcv;T  ^  ,  where  q  is  the  speed  of  the  flow  and 
F(av,bv;cv;r)  is  the  hypergeometric  function*  On  the  other  hand,  the 

particular  integral  of  the  velocity  potential  will  have  only  an  extra 
factor  involving  the  logarithmic  derivative  of  av,bv;cv;r^  •  Further¬ 
more,  "because  of  the  two  partial  differential  equations  connecting  them, 
"both  functions  will  have  common  coefficients.  Thus,  only  one  group  of 
hypergeometric  functions  and  one  set  of  unknown  constants  are  necessary 
for  the  complete  determination  of  the  two  functions. 


Construction  of  Solution  in  Holograph  Plane 

Before  outlining  the  procedure  adopted  here,  it  is  perhaps  proper  to 
describe  the  mapping  of  an  incompressible  flow.  In  order  to  simplify 
the  argument,  the  body  has  been  assumed  to  be  symmetrical  with  respect 
to  the  coordinate  axes  and  the  flow  at  infinity  to  be  parallel  to  the 
major  axis.  Such  a  flow,  when  mapped  onto  the  holograph  plane,  will  give 
rise  to  two  branches  of  Kiemann  surfaces,  each  corresponding  to  one— half 
of  the  physical  plane.  Because  of  symmetry,  the  discussion  may  be 
restricted  to,  say,  the  left— half  plane  D  (fig*  1)*  Because  of  the  fact 
that  the  maximum  and  minimum  velocities  occur  on  the  surface  of  the  body, 
the  whole  field  of  flow  is  mapped  into  the  interior  of  the  holograph  D 
the  boundary  of  which  corresponds  uniquely  to  the  lefi>-half  portion  of 
the  body  M’SM,  and  infinity  corresponds  to  a  point  P  on  the  positive 
axis  of  real  values.  Evidently,  this  point  is  singular. 

With  such  a  domain  in  the  holograph  plane,  it  is  possible  to  construct 
a  compressible  flow  from  a  "similar"  incompressible  flow.  By  similar  flow 
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is  meant  that  at  the  limit  of  zero  Mach  number  the  compressible— flow 
solutions  will  reduce  uniquely  to  the  given  incompressible  flow.  First, 
since  the  stagnation  point,  that  is,  the  origin  S,  is  a  regular  point, 
there  exists,  a  Taylor  expansion  in  its  neighborhood.  As  the  singular 
point  which  corresponds  to  the  infinity  of  the  physical  plane  is  an 
interior  point,  it  must  he  on  the  circle  of  convergence  C  of  the  solution. 
Second,  the  solution  represented  hy  the  Taylor  expansion  should  be 
continued  analytically  for  the  whole  domain.  This  can  be  done,  of  course, 
either  formally  by  analytic  continuation  or  by  solving  Cauchy’s  initial 
value  problem.  The  former,  however,  is  not  practicable.  As  to  the  latter, 
use  can  be  made  of  the  known  character  of  the  singularity  of  the  incom¬ 
pressible  flow  to  determine  the  form  of  the  expansion  outside  the  circle 
of  convergence  C.  If  the  Taylor  expansion  is  regarded  as  given,  the 
"outside"  solution  can  be  uniquely  defined  by  the  conditions  of  continuity 
across  the  circle  of  convergence  C.  This,  of  course,  is  based  on  the 
assumption  that  the  character  of  the  singularity  is  unchanged  by 
compressibility. 

It  Is  found,  however,  that  the  assumption  is  not  valid.  Theoretically, 
it  can  be  shown  that,  if  a  Taylor  expansion  corresponding  to  incompressible 
flow  is  given,  then  after  each  term  of  the  expansion  is  multiplied  by  a 
proper  hyperge  cane  trie  function,  the  resulting  solution  will  have  a  logar- 
rithmic  singularity  in  addition  to  those  it  originally  possessed.  This 
means  that  the  "distortion"  due  to  compressibility  becomes  larger  for 
larger  speeds.  It  can  also  be  shown  that  the  character  of  the  singularity 
is  preserved  for  nonvanishing  Mach  number  if  and  only  if  the  ratio  of  the 
specific  heats  assumes  the  value  of  -1.  In  order  to  take  this  effect 
into  account,  the  procedure  is  either  to  start  with  a  special  body  to 
compensate  it  or  to  eliminate  it  by  adding  an  extra  tern  to  the  outside 
solution.  The  former  procedure  is  difficult.  By  following  the  latter, 
the  part  of  the  outside  solution  which,  in  the  limit,  reduces  to  that  of 
the  incompressible  flow  is  regarded  as  known.  Then  the  Taylor  expansion 
and  the  extra  term  added  to  the  outside  solution  can  be  determined  by  the 
conditions  of  continuity.  The  extra  term  becomes  zero  for  zero  Mach 
number.  Thus,  it  is  seen  that  the  Chaplygin  condition  is  again  satisfied 
but  the  procedure  haB  been  greatly  extended. 

The  flow  with  a  finite  circulation  has  also  been  considered.  In 
order  to  simplify  the  mathematical  problem,  the  circulation  has  been 
assumed  to  be  very  small.  Under  this  assumption  the  effect  of  circulation 
on  the  solution  of  the  incompressible  flow  can  be  represented  approximately 
by  an  arbitrary  combination  of  vortices  and  doublets  at  the  singularities 
of  the  holograph.  Therefore,  to  the  original  singularities  there  will  be 
superposed  a  logarithmic  one  in  order  to  define  a  circulation.  Since  the 
added  part  due  to  circulation  is  an  even  function,  the  resultant  solution 
will  be  unsyrametric  with  respect  to  the  major  axis. 

The  solution  for  the  similar  compressible  flow  can  be  constructed 
in  the  same  manner,  with  the  exception  that  the  conditions  of  continuity 
across  the  circle  of  convergence  are  insufficient  to  determine  all  the 
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arbitrary  constants,  in  particular  those  characterizing  the  strengths  of 
the  vortices  and  the  doublets.  These  are  determined  from  the  condition 
that  at  the  stagnation  point  dx  =  0  and  dy  =  0  and  from  the  geomet¬ 
rical  condition  of  symmetry.  The  former  condition  is  satisfied  in  all 
cases  irrespective  of  whether  there  is  circulation,  whereas  the  latter  is 
required  only  when  the  circulation  is  finite. 


Improvement  of  the  Convergence  of  Power  Series 

The  whole  problem  thus  hinges  on  the  method  of  carrying  out  the 
actual  computation.  It  should  be  pointed  out  that  inasmuch  as  the  solu¬ 
tions,  in  the  limit,  reduce  to  harmonic  functions,  the  convergence  of  the 
power  series,  especially  in  the  neighborhood  of  the  circle  of  convergence, 

is  generally  very  slow.  This  is  particularly  true  of  the  part 
00 

q  VF/ay  —  v,bv  —  v;l  —  V;^  sin  v0  of  the  outside  series,  because 
0  V  y 

it  contains  the  hypergeometric  functions  which  increase  rapidly  with  v. 

This  situation  can  be  eased  somewhat  by  introducing  the  asymptotic 
expansions  of  the  hype rgeame trie  functions.  For,  after  these  are  substituted 
in  the  solutions,  the  first-  and  sometimes  even  the  second— order  teiro s  can 
be  summed.  As  a  result,  the  solution  in  each  case  can  be  broken  up  into 
two  parts,  one  of  which  is  of  closed  fona  and  the  other  is  a  power  series 
with  Improved  convergence.  Owing  to  the  fact  that  the  dominant  terns  give 
excellent  approximation  in  the  domain  of  validity,  the  value  of  the  tern 
given  by  the  power  series  is  usually  of  inferior  order.  For  practical 
purposes  when  high  accuracy  is  not  desired,  the  amount  of  labor  involved 
can  be  greatly  reduced. 

As  was  pointed  out  in  reference  8,  the  summed  part  can  be  identified 
as  the  "speed  distortion"  in  the  subsonic  region;  in  the  siipersonic 
region  where  the  differential  equation  changes  its  type,  it  can  be  inter¬ 
preted  as  a  "standing  wave,"  depending  only  on  the  two  characteristic 
parameters.  In  this  case  both  these  simple  solutions  are  known  to  be 
inaccurate,  especially  in  the  neighborhood  of  sonic  speed.  For  full 
discussions,  see  reference  8. 

It  must  be  added  that  in  the  case  of  the  derivatives  of  the  hyper- 
geometric  functions  the  dominant  terms  cannot  give  as  good  an  approxi¬ 
mation,  as  can  be  seen  by  comparing  f(^)(T)  with  in  table  1. 

For  this  reason,  the  coordinate  functions  which  involve  the  hypergeo- 
metric  functions  as  well  as  their  derivatives  have  more  important  cor¬ 
rection  terms  and  hence  require  the  use  of  many  more  terms  for  the  actual 
computations. 
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II  -  CONSTRUCTION  OF  A  SOLUTION  IN  HODOGPAPH  PLANE 
Transformed  Differential  Equations  and  Their  Particular  Solutions 


Let  u  and  v  be  the  velocity  components  of  a  two-dimensional  flow, 
parallel  respectively  to  the  x—  and  the  y— axis  of  a  Cartesian  system.  In 
the  case  of  steady,  irrotational,  and  isentropic  motion  of  an  inviscid, 
nonconducting,  and  compressible  fluid,  the  Eulerian  equations  can  be  inte¬ 
grated  to  give  the  pressure  p,  the  density  p,  or  the  sonic  speed  c 
in  terns  of  the  flow  speed  q: 


where  pc,  p0,  and  c0  denote  respectively  the  values  of  p,  p,  and  c 

at  the  stagnation  point,  and  y  is  the  ratio  of  the  specific  heats  of 
the  gas.  Furthermore,  because  of  the  kinematic  conditions,  there  exist  a 
velocity  potential  cp  and  stream  function  defined  by 


Pu  Po^y 
pv  =  Po*x 


where  the  subscripts  indicate  the  partial  derivatives. 


8 


NACA  TN  No.  Ikk^ 


By  use  of  equations  (k)  and  (5),  the  partial  derivatives  in  the 
expressions  for  dq)  and  cty  can  he  eliminated  so  that  dcp  and  d+ 
can  he  expressed  linearly  in  terns  of  dx  and  dy  or  dx  and  dy  can 
he  expressed  linearly  in  tenns  of  dcp  and  d\[r.  Furthermore ,  if  the 


Jacobian  function 


is  finite  and  nonvanishing,  the  correspondence 


between  the  physical  xy-  and  the  hodograph  ub— plane  is  1  to  1.  Under 
this  condition,  by  regarding  u  and  v  as  independent  variables,  the 
relations  connecting  the  differentials  yield  the  following  equalities: 


and 


xe 

ye 


e<pe  ■ 
e<Pe  + 


Po 


sin  6ife 


y  cos  ete  ) 


6<Pq 

®9q 


£2  sin 
P 


Po 

+  —  cos 


(6) 

(7) 

(8) 

(9) 


where  0  is  the  inclination  of  the  velocity  vector  to  the  x-axis  and 
accordingly 


u  =  q  cos  0 
v  =  q  sin  0 

Finally,  the  condition  of  integrability  demands  that 

Wq.  =  -  y  i1  -  m2)*0  U°) 

<Pe  =  — -c^,  (11) 

p  * 

where  M  =  q/c  is  the  local  Mach  number.  These  fundamental  systems 
permit  the  complete  determination  of  the  functions  <p  and  t.  For 
by  eliminating,  say,  <p(q.,0)>  the  resultant  .equation  for  +(q,0)  ia 
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when  v  is  a  positive  integer «,  where 


10 
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The  hype  rge  erne  trie  equation  thus  possesses  two  distinct  families  of  fundamental 
solutions  according  to  whether  v  differs  from  or  is  equal  to  an  integer.  The  second 
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integral  defined  by  expression  (l6)  was  found  most  appropriate  for  these 
particular  parameters  an  and  bn.  The  reason  for  subtracting 

Cn[+(bn)  —  Kcn)]  ■^(an>^n^cn’T)  from.  the  second  integral  defined  in 

reference  8  is  to  neutralize  the  contribution  of  the  first  integral  and 
the  discontinuities  carried  by  +(bn),  as  bn  is  negative. 

In  the  following  discussions,  the  two  fundamental  solutions  of 
equation  (13)  will  be  denoted  by  qvFv  (t  )  and  q“vF_v(T)  when  v  is 
not  an  integer  and  by  t)  and  q-^F_n(r)  when  V  is  an  integer, 

where  F_n(r)  is  defined  by  the  expression  within  the  bracket  in 

expression  (l6).  The  normalization  has  been  chosen  for  a  continuous  passage 
of  a  compressible  to  an  incompressible  flow.  The  most  important  property 
of  F_n(T }  from  expression  (16)  is  that  when  n  =  1 


F_x(t)  *1 

For  this  parameter  the  first  integral  reduces  to 


(22) 


FpM  s  F(1,-Pj2jt)  = 


(P  +  1)t 


1  -  (1  -  t) 


P+1 


(23) 


The  particular  integrals  of  equation  (12)  are  then  given  by 


<l%(T)jsinive 


when  v  is  not  an  integer  and  by 


(24) 


(25) 


when  v  is  an  integer  n.  In  addition  to  these  solutions ,  there  are 
two  other  integrals,  each  of  which  is  a  function  of  only  one  variable. 
On  assuming  if  =  ijr(q)  or  ty(e),  equation  (12)  yields 


(1  -r)P  £4 

T 


6  and 


j 


(26) 
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Since  the  particular  integrals  of  i|r(c,0)  are  known,  those  of 
<p(q,0)  are  shown  to  be  (see  appendix  A) 


(1-  T)-VFv(T)iv(T)j^Jve 

(1-  T)-Pq.-VP_V(T;  tv(r)  jflsj  V© 
when  v  is  not  an  integer  and 


(i-  TrVrn<T)!n(T)j^jBe 


(27) 


(28) 


when  v  is  an  integer  n.  Here  the  functions  lv(T)  and  ?_n(T) 
for  any  positive  value  of  v  are  defined  by 


V|v(r )  =  2t  logg  tV^v(t) 
v|_v(t)  =  2t  a.  ioge  t_V/2F_^(t) 


(29) 


The  following  expressions  corresponding  to  expression  (26)  are  similarly 
found: 


(1  -  T)-*  _  | 


(1  -  t)' 


dr 

T 


and  6 


(30) 


Hydrodynamic  Functions  of  Incompressible  Flows 

By  following  the  procedure  adopted  in  reference  8,  the  analysis 
starts  with  the  functions  required  in  defining  an  irrotational  incompressible 
flow.  In  the  case  of  an  inccmpressible  fluid  for  which  the  sonic  speed  is 
infinite,  the  equations  satisfied  by  the  velocity  potential  cp  and  the 
stream  function  i|r  become  haimonic.  If  W0(z0)  is  the  complex  potential, 
it  can  be  shown  that 
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Wo(z0)  =  <P0(xo^o)  +  ^o^o^o)  wi'fck  zD  =  x0  +  iy0  (31) 

In  a  simply  connected  domain,  the  functions  cp0  and  *  are  single¬ 
valued  and  continuous;  (pD  can  he  many-valued  only  when  the  flow  is 
now  free  from  circulation  rQ. 

If  w  denotes  the  ccmplex  velocity  u  -  iv,  it  connects  with 

Vzo)  by 

v  =  "¥Zq)  =  v(zo)  (32) 

dz0 

This  establishes  the  relation  between  the  physical  and  the  hodograph 
plane.  The  inverse  transfoimation 


zo  =  z0(w)  (33) 

exists,  provided  that  w*(z0)  ^  Q,  This  function  plays  an  important  role 
in  the  present  scheme  of  solution  and  will  be  known  as  the  transition 
function.  In  general,  as  it  is  an  inverse  function,  it  is  not  single- 
valued,  as  discussed  in  reference  8.  By  introducing  this  relation  into 
equation  (31),  the  complex  potential  in  the  hodograph  plane  is 


wo(w)  =  <P0(u*y)  +  i+0(u>y)  (3M 

In  case  the  solution  of  equation  (32)  is  many-valued,  WQ(w)  will 

represent  only  one  branch  of  its  many  solutions.  When  the  flow  is  not 
free  from  circulation,  the  function  w(zQ)  will  contain  r0  as  a 

parameter.  This  extra  tern  generally  makes  the  transition  function 
z0(w)  more  complex  than  it  would  be  if  rQ  vanishes.  In  practice,  the 

complication  can  be  reduced  to  a  certain  extent  by  a  linear  superposition 
of  the  two  effects  such  that  the  transition  function  becomes 


zo  =  zo(°)(v)  +  f2Zo(l)(w) 

4* 


(35) 


This  simplification  can  be  justified  as  long  as  the  circulation  is  weak 
so  that  teims  of  higher  order  may  be  neglected.  With  the  transition 
function  so  defined,  the  complex  potential  may  be,  similarly. 
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Here  zQ(0)(w)  and  WQ(0)(w)  are  respectively  the  transition  function 

and  complex  potential  for  zero  circulation;  z0(^)(v)  and  W0(^)(w), 
which  represent  the  effect  of  circulation,  are  known  when  WQ(w)  is 
given.  It  can  be  shown  that  generally  W^^w)  can  be  represented  by 
vortices  and  doublets  at  the  singularities  of  zQ(w).  To  this  order  of 
approximation,  it  is  easy  to  see  that  the  circulation  T0  is  correctly 
defined. 

Conversely,  when  W0(w)  is  given,  zQ(v)  can  be  obtained  by 
integration: 

z0(w)  =  J  W0' (w)  ^  +  Constant  (37) 

which,  in  fact,  is  equivalent  to  equation  (32). 


Construction  of  a  Symmetric  Solution  about  the  Origin 

From  the  considerations  of  reference  8,  if  the  flow  about  a 
symmetric  body  placed  at  the  origin  of  the  xy— plane  is  mapped  onto  the 
hodograph  plane,  the  whole  left-half  plane,  exterior  to  the  body,  will 
correspond  to  a  region  in  the  hodograph  plane,  of  which  the  point 
w  =  U,  U  being  the  modulus  of  w  at  infinity  of  the  z0-plane,  is 
a  singularity.  Then  the  domain  within  the  circle  |w|  =  U  is  single- 
valued  and  regular.  If  the  complex  potential  WQ(w)  is  associated  with 
a  definite  flow  in  the  zQ-plane,  it  must  be  analytic  and  regular  within 
|w|  =  U.  Consequently,  it  has  the  Taylor  expansion: 

00 

Wo(w)=-^Anwn  |w|  <  U  (38) 

0 

where  the  coefficients  An  are,  in  general,  complex.  If  the  body 

is  symmetrical  with  respect  to  both  coordinate  axes,  then  the  coefficients 
are  real.  Separating  into  real  and  imaginary  parts  yields,  according  to 
equation  (3% 


00 


+  0(<l>e)  =  3>~  Vin  Sin  n0 

(39) 

2 

00 

q  <  1 

<PO(<1,0)  =  &n<Ln  C0E 

(40) 
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where  w  =  qe“ i0  and  A^  =  0  because  w’(z0)  ^  0  at  w  =  0.  Frcsn 
now  on,  zQ,  w,  p,  and  p  are  normalized  in  terms  of  a,  U,  p0,  and 

p0,  respectively.  Then  q  =  1  at  infinity,  and  p  and  p  will  be 
unity  at  w  =  0. 

According  to  Chaplygin’s  procedure,  the  corresponding  solutions 
for  the  compressible  fluid  can  be  obtained  by  simply  replacing  the 

function  qn  in  equations  (39)  and  (40)  by  qnFn(r)(T)  and  - 

(l  —  T)”PqnFn(r)(T)?n(T),  respectively,  as  shown  by  expressions  (25)  and  (28) 

The  second  integrals  are  excluded  by  the  condition  of  regularity  at 
q  =  0.  Thus  the  following  equations  are  obtained: 


=3EI  sin  n9  (41) 

2 

q  <  1 

00 

cp(<l,e)  =  -  (1  -  An^n^(T)^n(T)  008  ne  +  Constant  (^2) 

2 


vhere 


F  (p)(T)  = 

n  Fn(Tl)  •^(an,bnicn»T1) 


m 


and 


T 


c  1  u2 

1  2P  CqS' 


that  is. 


t 2_  corresponds  to  the  free-stream  velocity  U. 


It  is  seen  that  if  cQ  ®,  then  both  t  and  T]_  tend  to  zero  and 
Fn(r)(r)  -»  1.  Thus,  the  solutions  are  reduced  to  the  incompressible  form. 

Furthermore,  if  q  1,  the  character  of  the  solution  is  exactly  like 
that  of  the  incompressible  solution.  Hence  all  the  specified  conditions 
are  satisfied.  It  must  be  remembered  that  these  conditions  are  valid  only 
for  subsonic  flows,  and  for  this  reason  ti  is  restricted  to  the  subsonic 
region. 


The  series  (equation  (4l))  constructed  in  this  manner  is  actually 
convergent  and  represents  the  function’  \jr(q,0)  within  the  circle  of 
convergence  q  =  1  (reference  8). 


Analytic  Continuation  of  Solution  (Branch  Point  of  Order  1) 

In  this  section,  it  is  proposed  to  continue  the  solutions  t  and  cp, 
represented  respectively  by  equations  (41)  and  (42),  analytically  outside 
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the  domain  |w|  <1.  The  domain  outside  |w|  <  1  is  generally  many- 
valued.  In  order  to  be  precise,  let  it  be  a  branch  point  of  order  1. 
Generally,  the  function  W0(w)  has  other  singularities  in  addition  to 

the  one  at  w  =  1.  However,  such  singularities  lie  outside  the  domain 
of  interest  and  thus  need  not  be  investigated.  Let  the  nearest  singularity 
be  given  by  w  =  V  >  1.  Then,  the  domain  to  be  considered  outside 
|w|  =  1  is  an  annulus  with  a  cut  joining  the  two  singularities.  The 
proper  representation  of  WQ(w)  in  such  a  region  which  has  a  branch 
point  of  order  1  at  w  =  1  is 


_  1 

WQ(v)  =  iw  2  WQ*(w)  (44) 


where  WQ*(w)  is  single-valued  and  regular  within  the  open  annulus 

1  <  |  w  |  <  V.  Hence,  in  any  closed  domain  1  +  5<  |w|  <  V  —  &,  &  being 
a  positive  value,  there  exists  a  uni f orally  and  absolutely  convergent 
series : 

W0*(w)  =  X  (Vw*  +  Cn*w-n)  (45) 

which,  on  substituting  in  equation  (44),  will  give  the  continuation  of  the 
Taylor  series  (equation  (38)).  This  is 

W0(w)  =  i  21  ( Bnwv  +  Cnw“v  ^  1  <  | w |  <  V  (46) 

where  the  constants  Bn  and  Cn  are  real  because  of  symmetry  of  the 
body  and  v  =  n  +  rj,  n  being  a  positive  integer. 

In  order  to  continue  the  solution  of  +(<i>e)  outside  the  circle 
g.  =  1,  two  alternatives  are  encountered*  Suppose.,  first,  that  the 
character  of  the  singularities  of  t|ro(q,0)  defined  by  equation  (39)  iB 
•unmodified  by  the  hypergeometric  functions*  Then  the  solution  for  the 
compressible  fluid,  valid  in  the  annulus  1  <  |w|  <  V,  can  be  obtained 
by  introducing  the  proper  hype rgeorne trie  functions  corresponding  to  the 
parameter  v.  The  continued  solution  would  be 

M<1>0)  =  X  [®nTV?v^r^(T)  +  Cnq_vF_v(r^(r)  j  cos  v0  1<  q<  V  (47) 


where 
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I\Mt) 

F.„(r,(T)  = 


_  Fy(r) 
~ 

F— y(T) 
F_v(ti) 


> 


(48) 


Here  FV(T)  and.  q  ^^(t)  are  respectively  the  first  and.  second. 

A/ 

integrals  of  the  hypergecmetric  equation;  Bn  and.  Cn  are  constants. 
It  may  he  added,  that  the  coefficients  Bn  and  Cn  are  not  the  same  aB 

those  in  equation  (46)  for  the  incompressible  flow  but  should  be 
determined  by  the  conditions  of  continuity. 

Since  the  partial  differential  equation  considered  here  is  of  the 
second  order,  in  order  to  ensure  that  the  analytic 

continuation  of  1rin(q,0),  two  conditions  must  be  satisfied  at  the 

boundary  of  the  respective  regions  of  convergence,  that  is,  the  circle 
q  =  1.  The  two  conditions  are  as  follows: 


*in^0) 

q=l 


=  W*'0) 


q=l 


(49) 


^in(<l,0)  q=1  =  £5* out(^0)  q=i 


(50) 


It  should  be  noted  that  the  condition  (equation  (49))  in  this  particular 

Vin(*e)|  d 


case  is  identical  with  |^4rln(q,0)  ^  =  |g'tfOut(^-’0)  J  _1  and  that  the 


values  taken  on  the  circle  of  convergence  are  the  limiting  values  if  such 
limit  exists  in  each  case^  These  conditions  would  be  sufficient  to 
determine  both  and  Cn.  However,  the  solution  constructed  by  this 

method  suffers  serious  distortion  because,  even  though  the  values, 
say,  \|r  and  ty0,  agree  on  the  circle  q  =  1,  the  tangent  of  4r  =  Constant 
on  q  =  1  deviates  from  that  of  \|r0  =  Constant.  This,  of  course,  will 
affect  the  character  of  the  function  tcXl-*0)*  If  W0(w)  defined  by 

equation  (38),  the  question  arises  as  to  the  singularity  after  the  coeffi¬ 
cients  are  multiplied  by  the  proper  hypergeometric  functions.  The  exact 
answer  to  this  question  has  not  been  attempted,  but  a  rough  estimation 
given  under  theorem  1  seems  sufficient  for  the  present  discussion. 


Theorem  1.—  A  given  Taylor  expansion  such  as  equation  (38),  which 
has  a  singularity  at  w  =  1,  iB  modified  by  multiplying  its  coefficients 
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~  1 
n^n  = 


. 

r 

ri 

hI„(ti)  -  n|  (ti) 

\° 

L 

n  !n(ri)  -  i_n(Ti)  ^  0 

h5_^(ti)  -  n|_n(ri)  |  j>  -  2jtD0i_n(rl) 


®m 


(5D 


„  .  £  -j  I  -  n|  (ti) 

ln(n)  -  l_n(Tl)  \  0  1  I  ^ 


nDn  -i - i 


•*'11  -  n?n(Tl> 


°m  -  2lSo*n<Tl> 


(55) 


where  the  determinant  n(£n  -  i_n)  ^  0.  It  is  observed  that,  as 

fs* 

t-,  0,  Dn  -»  0  and  An  An.  Consequently.,  Chaplygin’s  condition  is 

again  satisfied.  Furthermore,  1^=0  also  for  7  =  —1.  In  that 
event,  the  compressible  and  incompressible  f lows  will  have  the  same 
singularities. 

The  solution  is  formal.  In  order  to  prove  that  the  function  f(q,0) 
represented  by  equations  (ll)  and  (51)  is  regular  in  the  respective 
domains  of  validity,  the  truth  of  the  following  theorem  must  first  be 
demonstrated. 

Theorem  2.-  If  the  constants  An  and  are  defined  by  equations  (5!) 

and  (55),  respectively,  and  if  the  series  (equations  (38)  and  (16)) 
converge  uniformly  and  absolutely  in  the  danains  specified,  the  series 
(equation  (11))  with  coefficients  An  and  equation  ( 51)  are  uniformly 
and  absolutely  convergent  in  the  corresponding  domains  (appendix  E). 

With  ^(q,9)  so  determined,  the  velocity  potential,  as  given  in  the 
section  entitled  "Transformed  Differential  Equations  and  Their  Particular 
Solutions,"  is 


00 

cp(q,e)  =  -(1  —  r)-^  AnqnFn(r)(T  }  cos  nO  +  Constant  (56) 

2 


and 
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repre  s  entati on : 
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This  assumption  is  Justified  only  by  the  fact  that  ro  is  small  and  the  corresponding 
cp^1^(q.,0)  yields  a  correct  circulation.  In  the  case  of  a  symmetrical  body,  t^°^(q.,0)  was 
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Here  the  constants  A^l),  Bj^-1-),  and  Dn(^  characterize  the 

strength  of  the  vortices  and  doublets  and  are  determined  by  the  conditions 
of  continuity  at  q  =  1.  Generally  P(q,0)  is  a  polynomial  and  is 

introduced  to  add  to  i|r^(q,0)  a  mixed  symmetry.  The  number  of  terms 
required  is,  from  cause  to  effect,  unknown  and  may  be  different  for 
different  bodies.  In  the  present  case,  only  one  term  is  taken  for 
simplicity.  On  assuming  that  every  series  converges  uniformly  and- 
absolutely  in  the  respective  domain  of  validity  and  tends  to  a  definite 
limit  as  q  1,  the  conditions  of  continuity  at  q  =  1  are 


+  cn(1)j’_n(T1) 


cos  n©  =  —A, 


H  n 
1 


AnC1)Fn(r1)|n(Ti)  +  CnU)F_n(TiH_n(Ti) 


(1) 

I 

o  <  e  <  2* 

cos  nB  =  0 


and 


£  [Bn(l)Fn(ri)  _  D^^F.Jti) 


cos  n©  =  0 


0  <  ©  <  2jt 


II  n  [Bn(1)Fn(r1)ln(T1)  _  Dn(l)F_n(Ti)l_^(xi)  j  cos  n©  =  -®oU)(l 


By  the  uniqueness  theorem  of  expansion,  since 


cos  n©  =  — 0  <  0  <  2n 
1  d 


it  follows  that 


An(1)Fn(rl)  +  cJ^F^ti)  =  2A0(1) 


An(1)Fn(T1)|n(T1)  +  cJ1)F_n(T1)|_n(r1)  =  o 


and 
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Bn(l)Fn(Tl>  “  Dn(l)F-JTl)  =  0 

Bn^Fn(Tl)Sn(Tl)  “  Dn^F_n(TlH-n(Tl)  =  f  U  -  Tl)^1* 

The  solutions  are,  respectively: 

An(l)  =  -A0(l)(l  "  Tl)"^F-n(TlH-n(Tl) 

Cn(l)  =  Ao(1)(l  -  Tir^n(TlHn(Tl) 

B  (1)  -  B  (1)  *=o££ 


n 


(1)  s  b  (1)  Mu 1 


Vx;  =  Bo 


(67) 


(68) 


From  theorem  2,  it  can  easily  be  shown  that  'i'(^)(q.,0)  so  defined,  is, 
as  assumed,  uniformly  and  absolutely  convergent  in  any  closed  domain  in 

q  <  1  and  1  <  q  <  Y  and  that  =  0  when  q  =  0.  Furthermore, 

the  arbitrary  constants  A0(^-),  B0(^)>  and  CQ^)  are  to  be  determined 

by  the  auxiliary  conditions.  In  the  first  place,  the  condition  that 
dx  =  dy  =  0  at  q  -  0  demands  that 


-A^1)  +  B^1)  =  0 


which  leads,  because  of  expression  (22),  to 

A0(1)  =  (1  -  t1)Pb0(1)  (69) 

In  case  there  is  no  stagnation  point,  such  as  in  the  case  of  a  body  with 
a  sharp  trailing  edge,  this  condition  is  again  to  be  satisfied  in  order 
to  avoid  the  multiplicities  in  x  and  y.  The  remaining  two  constants, 
namely,  BqW  and  CQ(1),  will  be  determined  by  the  condition  of 
symmetry  and.  will  be  discussed  under  Determination  of  Integration  Constants. 
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Similarly,  the  function  q>(q,0)  is 


<P(l,9)  =  <P(o)(4,0)  +  TAp(1)(ct.e) 

4it 


(70) 


Here  the  function  <p(°)(<lj0)  is  given  hy  equations  (56)  and  (57),  and 
^■^(cbO)*  for  q  <  1,  is 


<p(l)(4,8)  »  <1  -  r)-S 


JZ  4i*1*l”Fn(T)!n(T)  sln  ”6  +  Q(a,8) 

2 


+A0  (71) 


and,  for  1  <  q  <  V,  is 


<p(1)(q.,0)  -  (1  -  T)-*' 


®o(l)  Z  Bn(1)q1¥n(T)|11(T) 
1 


+  Cn(l)q-nF_n(T)^_n 


sin  n0  +  Q 


l  -  B0(1)(n  +  0) 


(72) 


vhere 


Q  =  -C0(1)q.?P2(T)S2(T)  cos  20 

*n(1)  -  *o(1)  [  U<T1>  +  s]  F-*<T1>  +  Ao'1^1  -  J) 

bJ1)  =  (1  -  Tl)P(^l  -  £) - - - 

V  n/Fn(Tl)Tn 

°n(1)  =  [^(n)  +i]ra(Tl) 

The  velocity  potential  thus  gives,  hy  definition,  a  circulation  T: 

r  =  W  dq>  =  r  cpj1) 


r  (73) 


(7*0 


where  the  integral  is  taken,  in  clockwise  direction,  about  the  branch 
point  in  two  complete  circuits. 


When  a  pair  of  functions  f(q,0)  and.  <p(q,0)  are  constructed  in  the  hodograph  plane 
the  flow  pattern  to  which  they  correspond  in  the  physical  plane  cannot  "be  given  directly 
hut  further  integration  of  two  differential  equations  is  required.  This  process,  however, 
does  not  involve  any  difficulty  owing  to  the  fact  that  both  dx  and  dy  are  total 
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Differentiating  equations  (75)  (76)  partially  with  respect  to  q  and  using  equation  (13) 

to  eliminate  the  derivative  higher  than  the  first  yields  x^  and  y^,  By  comparing  these  tw 

expressions  with  those  obtained  directly  frcm  equations  (8)  and  (9)  it  can  be  concluded  that 
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In  the  case  of  circulatory  flow  corresponding  to  equations  (63)  and  (70),  the 
transition  functions  x(q,0)  and  y(q,0)  can.  he  written  as 


30 


NACA  TN  No.  1445 
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cos  Q  +  yout 
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Here,  since  Bq(^)  =  0,  the  terns  for  n  =  1  again  contribute  no  diffi¬ 
culty.  The  constants  of  integration  in  equations  (83)  and  (87)  are  left 
out  because  they  can  be  incorporated  in  equations  (75)  and  (77) > 
respectively. 


Determination  of  Integration  Constants 

Consider  first  the  case  when  the  flow  is  symmetrical.  The  transition 
functions  x(q,0)  and  y(q,0)  for  such  a  flow  are  given  respectively  by 
equations  (75)  and-  (76)  for  q  <  1  and  by  equations  (77)  and  (78)  for 
1  <  q  <  V,  involving  four  arbitrary  constants.  Of  these  four  constants, 
xout  and  y^n  can  be  chosen  arbitrarily  by  translating  the  coordinate 

axes.  Indeed,  because  of  symmetry,  it  is  preferable  to  choose 


xout  (<1*0)  =  0  for  <1  =  <3u  811(1  0=0  (89) 

and 

yin  (<1>0)  =  0  for  q  =  0  (90) 


whence 


'out 


-  ~-Th>-  ■■  J  E 
%  1  n  "  ^  L 


(91) 


yin  =  0 


(92) 


where  qy  denotes  the  flow  speed  on  the  upper  surface  at  the  midsection  of 
the  body.  Then  xou^  depends  on  and  tends  to  zero  with  t^. 

The  other  two  constants  can  now  be  determined  by  the  condition  of 
continuity  at  q  =  1,  namely. 


xin(i,e)  =  xout(i,e) 


(93) 

(94) 


yin(!,e)  =  youtf1*0) 


0  <  0  <  2n 


6 


HACA  TIT  No.  1445 


33 


Since  both  limits  exist,  the  integration  of  the  identities  yields 


*in 


=  jl  ZJlL  t  J 

2n 


+  H  v(Bn  +  cn) 
0 


II  V  [sn^Ti)  +  Cn^T!) 
f__i _ i_ 

|_(v-l)2  (v  +  1)2 


1  +  1 


(v  -  l)2  (v  +  1) 


+  X 


out 


y out  =  ® 


(95) 


(96) 


The  fact  that  yout  =  0  is  the  consequence  of  y(q,0)  =  -y(q,2it  —  0); 
2xj_n  will  be  defined  as  the  chord  of  the  body. 

For  the  flow  with  circulation,  according  to  the  preceding  section, 
there  Eire  six  constants  instead  of  four.  Let  those  four  arising  from 
integration  be  considered  first.  Corresponding  to  equations  (89)  and  (90) 
the  following  equations  may  be  chosen: 


■out 


(o)  _  (j  -TU)  P 


<3U 


II  ~2 - nF-n(Tu)  [n^-n(Tu)  +  1 

1  -  1  L  j 


+  50  “  ^px(^I»°)  f 


7  out 


(1)  = 


(97) 

(98) 


Then  the  condition  of  continuity  at  q  =  1  gives 

r  v[Bnlv(Tl)  +  Cni^tr-L)] 


in 


(o)  =  Li- n rt  J 

2jt 


1  +  1 


(V  -  l)2  (v  +  l)2 


>  +  Iout(°) 


(99) 


(100) 
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In  order  to  derive  the  second  term  in  equation  (100),  use  has  been  made 

of  the  fact  that  F^(t)  -  =  (1  —  t)^,  which  is  just  the 

Wronskian  of  the  particular  integrals  of  the  hypergecmetric  equation 
for  n  =  1. 

The  arbitrary  constants  B^1)  and  CQ W ,  on  the  other  hand, 
are  determined  by  an  entirely  different  consideration.  The  fact  that 
bJ1)  and  C0^  are  different  from  unity  and  zero,  respectively,  as 
they  would  be  if  Tq  tends  to  zero  is  due  to  the  fact  that  the  distortion 

produced  by  compressibility  is  nonuniform  at  the  surfaces  of  the  body. 

In  order  to  correct  this  defect  completely,  a  more  elaborate  method 
would  have  been  required.  For  the  present  simple  investigation  in  which 
symmetry  is  not  strictly  satisfied  even  in  the  limiting  case  of  zero 
Mach  number,  the  condition  of  symmetry  will  be  applied  to  only  a  few 
selected  points.  First,  let  it  be  required  that 


Xout^'0)  ~  0 
Xout(V2Tr)  =  0 


(101) 


where  q^  stands  for  the  flow  speed  on  the  lower  surface  at  the  mid- 
section  of  the  body.  These  two  equations  detezmine  uniquely 
namely., 


£*  (i)  =  _ 

°  e2(qu)  -  e2(qL) 


(102) 


where 


sc,)  -  ^---T.ng 

<1 


H  -T^nd-^-nMC^-n  +  1)  +  5o 
1  nd  -  1 


(103) 


e2(q)  =  (1  -  t)_13<1F2(t)  [2?2(t)  +  l]  (104) 

The  value  of  CQ(^)  becomes  zero  with  rq  for  Dn  0  as  rq  -»  0,  as 
defined  in  equation  (55)* 

The  maximum  distortion  with  respect  to  the  y— axis  will  occur  at  the 
midsection.  If  the  required  body  is  assumed  to  be  symmetrical,  that  is. 


y(qy^O)  =  y(qL,2*) 
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the  shape  of  the  body.  If,  however,  the  camber  is  not  controlled,  A^ L>  may  be  taken  equal 

to  unity,  by  analogy  with  the  case  of  incompressible  flow.  Then  the  condition  of  conformality 
(equation  (69))  becomes  the  only  condition  to  define  the  circulation.  The  relative  merit  of 
these  two  procedures  can  be  ascertained  only  by  comparing  the  numerical  results. 
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III  -  IMPBOVEMENT  OF  CONVERGENCE  OF  SOLUTION  BY  ASYMPTOTIC 
PROPERTIES  OF  HYPERGEOMETPJC  FUNCTIOIB 
Trans foimat ion  of  Stream  Function 


The  stream  function  i|r(q,e)  for  a  flew  which  was  derived  from 
\lro(q.,0)  with  a  branch  point  of  order  1  has  been  given.  The  form  of 
representation  is  not,  in  general*  suitable  for  practical  calculation. 

The  difficulty  is  twofold:  First,  the  series  involves  an  infinite  number 
of  hypergecmetric  functions  which  are,  in  turn,  defined  as  infinite 
series.  The  convergence  of  the  hypergeametric  series  in  this  particular 
instance  decreases  with  increase  of  the  parameter  v.  This  means  that  the 
computation  for  the  later  terms  of  the  series  for  1f(q.,0)  will  be 
increasingly  laborious.  Second,  the  convergence  of  the  power  series 
defining  the  function  +(q,e)  is,  as  expected,  very  slow  in  the  neighbor¬ 
hood  of  the  circle  of  convergence.  In  order  to  render  the  method  of  any 
practical  value,  the  task  of  transforming  the  series  into  a  more  rapidly 
convergent  one  is  encountered.  For  this  purpose  the  following  procedure 
is  adopted. 

The  stream  function  +0(q,e)  for  the  similar  incompressible  flow 
is  (see  section  entitled  "Construction  of  a  Symmetric  Solution  about  the 
Origin" ) 


t0(^e)  =  £  Vi1*  sin  ne  1  <  1 

2 

which  is  absolutely  and  uniformly  convergent  in  any  closed  domain  in 
q  <  1.  Then  it  can  Justly  be  thought  of  as  representing  not  merely  a 
regular  but  a  closed  function.  In  doing  this,  of  course,  a  large  class 
of  problems  is  automatically  eliminated  and  those  cases  are  then 
considered  in  which  simple  representation  of  both 

exists.  This  is  Justified  only  by  the  mounting  difficulties  faced  in 
carrying  out  such  a  detailed  investigation. 

It  is  thus  observed  that  the  difference  between  the  stream  functions 
'Kq^e)  and  to(q,0)  lies  only  in  the  appearance,  in  the  former,  of  the 

hypergeometric  functions.  If,  however,  approximation  is  allowed,  then, 
according  to  expression  (223),  first  let  Fh(t)  be  substituted  by  its 

asymptotic  expression,  namely, 

*n(r,<u  ~  [ 1  +  °@] 


n  >  N  (107) 
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where  t(-r)  =  .  Here  the  tenn  involving  f(^)  (t  )  «  1  has  been 

T(t!) 

neglected.  Furthermore,  as  will  be  shown  in  equation  (159),  to  the  same 
order  of  approximation  the  coefficient  An  can  be  written  as 

*tl  ~  ^  +  °(h)  n  >  N  (l08) 

By  substituting  the  approximate  values  in  equation  (41)  it  can  be  shown 
that 


♦(4,0)  ~  +o(«b,0)  (109) 

f(Tl) 

That  is,  to  this  order  of  approximation  the  power  series  representing 

the  stream  function  +(q,0)  can  he  summed  and  is  given  by  expression  (109). 

As  was  shown  in  reference  8,  the  asymptotic  representation  is  valid 
only  when  the  parameter  n  is  large.  Namely,  during  the  summation  of 
expression  (109)  the  value  neglected  becomes  smaller  as  n  increases  and 
*  approaches  zero  as  n  tends  to  infinity.  This  concentration  of  errors 

in  the  lower-order  terms  makes  it  especially  easy  to  apply  the  correction 
if  high  accuracy  is  required.  In  doing  this,  the  quantity  given  by 
r-  expression  (109)  can  be  added  and,  at  the  same  time,  subtracted  from 

^(q,6).  Then  a  simple  manipulation  shows  that 

♦  (<1,6)  =  ^(q,©)  +  +2(4,0)  (HO) 

where 

■♦1(4,0)  =  ♦0(<it>e)  (m) 


♦2(4,0)  =  YL  Gn(T)4n  sin  n0 
2 


q  <  1  (112) 


with 


Gn(T)  -  Fn(T)AAn  + 


AFn(r) 


(113) 


BJa _ ia. 


Fn(Tl)  ^UJ^Ti) 

AFn(T)  =  Fn(r)  -  f(T)Tn(r) 


(114) 
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Here  n  is  a  positive  integer.  The  stream  function  t(q,0)  is  then 
represented  "by  the  sum  of  two  functions  +p(q,0)  and  Of 

these,  +1(11,0)  is  of  closed  form,  which  differs  from  +D(q i,0)  only 
by  a  change  of  scale  of  q,  and  +2(1^)  is  a  difference  of  two  absolutely 
and  uniformly  convergent  series  and  hence  is  absolutely  and  uniformly 
convergent.  In  fact,  according  to  expressions  (107)  and  (108),  G^(t) 

is  of  order  W>(  t)j  the  convergence  of  +(q,0)  is  therefore  increased 
by  l/n.  This  actually  is  the  gist  of  the  whole  problem. 

In  the  annulus  region  1  <  q  <  V;  on  the  other  hand,  the  stream 
function  tyo(q,0)  is  represented  by 

00  /  v 

+o(t>0)  =  ^BnqV  +  Cnq"V  J  cos  V0 
and  the  stream  function  +(q,0)  of  the  compressible  flow  is 

+(q,@)  =  J 
0 

+  Dc(*  -  ©)  +  f:  5n«l"nF-n(r)(T)  sin  n0 
1 

where  Bn  and  Cn  are  assumed  to  be  given  by  to(q,0),  and  Fv^Ct) 

stands  for  the  ratio  of  Fv(t)  to  Fv(ri).  In  this  region,  the  hyper- 

geometric  functions  will  be  of  mixed  character.  If  the  critical  point 

t„  =  - - —  is  not  reached,  they  are  of  exponential  type;  beyond  this 

c  2p  +  1 

point  they  will  change  over  into  oscillatory  type.  If  tc  lies  in  the 
range  1  <  q  <  Y,  the  singularity  of  the  asymptotic  expansions  of  the 
hypergeometric  functions  will  certainly  be  found  inside  the  domain  in 
question.  If  this  neighborhood  is  excluded  and  the  hypergeometric  functions 
are  first  substituted  by  their  asymptotic  forms  given  in  expressions  (223) 
and  (224),  it  can  similarly  be  shown  that,  for  T]_<  t  <  ., 

+(q>0)  =  +i(q>0)  +  +2^<i»e)  +  +3(<1j0)  (H5) 


BnqVFv(r)(T)  +  Cnq“VF_v(r)(r)]  cos  V0 


+l(<l,0)  =  +o(<l'0) 

i|r2(q,0)  =  ^  j^BnHv(T)qV  +  CnH_v(r)q“V  j  cos  v0  (ll6) 


where 
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*•3(1,0)  =  D0(it  -  0)  +  y~  Onl  nF_i/r^(T)  ain  "0  (117) 


and 


H^t)  =  Fv(t)AFv“1(t!)  +  - i- - AFv(r) 

f(Tl)T  (Ti) 

H_v(t)  =  F_v(t)AF_v_1(t1)+ - - AB1v(t) 

f(Ti)r-v(r1) 


>  (118) 


with 


£F  ”"1( t i)  =  — ± - - - 

V  Fv(ri)  f(Ti)Tv(Tl) 


AFv(t)  =  Fv(t)  -  f(r)TV(r) 


In  equation  (115)  +i(l,0)  again  represents  a  closed  function 
*o(lt»0),  and  *2(1,0),  8X1  absolutely  and  uniformly  convergent  series 

with  improved  convergence.  The  fact  that  ^3(1,6)  is  not  summed  is  due 

A/ 

to  the  fact  that  Dn  is  of  inferior  order  as  compared  with  An.  Even 
though  Dn  decreases  as  l/n,  *3(1,0)  will  require  fewer  terms  than 
*2(1,0)  although  its  coefficients  behave  like  if “3/2  (provided  that 

the  flow  *0(x0,y0)  Is  °1  the  nature  of  a  doublet).  Here  the  functions 
Hv(t)  are  so  defined  that  the  functions  of  Tp  and  those  of  t  are 
separated  so  as  to  make  possible  the  tabulation  of  AFv(T)  and  AF_v(t). 

In  the  region  - - —  <  t  <  1.  on  the  other  hand.  FV(T)  and 

2(3  +  1  v 

F_v(t)  are  replaced  respectively  by 

i  f(r)Tv(r)  cos  (vs)  -  •2) 

•|  f (t  )T”^  (t  )  cos  (vcd  +  £) 

where  f(T),  T(t ) ,  and  cd(t)  are  given  in  table  1.  The  factor  1/2 
is  introduced  before  the  first  expression  for  symmetry.  By  writing 


4o 
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2  cos  ve  cos  (vcc  -  ~  (cos  v|  +  cos  Vtj  +  sin  v!  -  sin 

2  cos  vO  cos  (vo>  +  A)  =  JL  (cos  v|  +  cos  vT)  -  sin  v|  +  sin  vt)  ) 

4  \[2  \  / 


with 


I  =  e  +  os 

T]  =0—03 


(119) 


the  following  expression,  corresponding  to  equation  (ill),  is  obtained 


♦(<1.>0) 


f(T) 


!/2f(Tl) 


♦o(^8)  +  ^(M)  +  <P0(MJ  “  <P0(M)] 


where  X  is  a  constant  defined  by 


X  = 


g(gp) 


a/ 2 


(1  +  a)a(2px1)1/2  T(tjl) 


>  1 


(120) 


as  from  eauation  (228)  qt  =  XU  if  -  <  t  <  1*  The  constant  X 

2  p  +  1 

then  is  a  function  of  the  free— stream  Mach  number  and  the  characteristic 
constant  of  the  gas  hut  independent  of  the  shar^e  of  the  body.  By 
eliminating  the  error  introduced  during  summation,  the  stream  function 
♦(q,0)  in  the  supersonic  range  —  -  <  t  <  1  is 


♦(q,0)  =  tx(q,e)  +  *2(q,0)  +  *3(q,0) 


where 

+l(4,0)  =  2-5/2  1 Jlly  3|f0(X,|)  +  i|r0(M)  +  <po(X,0  -  <P0(M) 


(121) 


♦2(4,0)  =  X!  f  BnH  (r)qv  +  CnH_v(r)q  v 
0  1 


cos  v6 


(122) 


Here  \|r3(q,0)  is  given  by  equation  (117)  and  ^(t)  and  H  (t),  by 
equation  (118),  except  that  AFv(T  )  and  £F_v(t  )  are  defined  by 
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£F  (t) 

±v 


=  F  (t)  _  If  (t)^  cos 
±v  2 


(123) 


Unlike  the  previous  calculations,  H  (t)  in  equation  (122)  is  not 
of  the  order  of  l/v,  owing  to  the  presence  of  l/2  in  front  of 

f (t )TV  cos  (va>  -  —)•  This,  however, '  does  not  introduce  a  serious 
objection,  as  the  series  in  which  it  appears  now  behaves  like  BnXv,  which, 
according  to  equation  (120)  (X  >  1),  converges  more  rapidly  than  BnqV 
in  *0(q,9). 

In  the  hyperbolic  domain,  moreover,  the  function  +  q(q,0)  depends, 

aside  from  a  factor  f(T),  only  on  the  two  independent  families  of  the 
characteristic  parameters  l  and  r\  defined  by  equation  (119).  This 
result  is  most  striking,  as  it  shows  that  the  main  part  of  the  solution 
satisfies  the  simple  wave  equation  and  thus  clearly  demonstrates  its 
hyperbolic  character.  With  both  the  incompressible  stream  function  to(q,0) 
and  the  incompressible  velocity  potential  cp0(^.^)  appearing  in  the 
solution,  it  is  impossible  to  establish  a  simple  relation  between  the 
incompressible  streamlines  and  the  compressible  streamlines.  Since  such 
a  simple  relation  is  the  foundation  of  the  so-called  speed  correction 
formula  for  a  qxiick  estimation  of  velocity  distribution  in  compressible  flew 
from  that  of  incompressible  flow  over  the  same  body,  this  idea  cannot  be 
extended  to  supersonic  regions.  On  the  other  hand,  this  also  indicates 
that,  although  the  differential  equation  for  \|/(q,0)  is  hyperbolic  in 
the  supersonic  range,  it  cannot  be  reduced  to  the  simple  wave  equation 
by  a  mere  distortion  of  the  speed  scale  as  given  by  the  function  cjd(t). 

For  if  this  were  the  case,  then  ^p(q,0)  would  constitute  an  exact  solution 
without  the  additional  ^(q^)*  This  fact  is  all  the  more  important  as 
the  additional  ^(q^)  is  not  small  in  comparison  with  typ(q,0)  for  the 

mixed  subsonic  and  supersonic  flows,  especially  for  the  transitional 
region  near  the  sonic  velocity.  However,  in  the  case  of  pure  supersonic 
flow,  +2(1^0)  might  be  small;  then  tp(q^)  alone  may  be  used  as  a 

satisfactory  approximation. 


Transformation  of  Coordinate  Functions  x(q,0)  and  y(q,@) 

From  equations  (75)  and  (76)  it  is  seen  that  the  coefficients  of  the 
series  defining  the  coordinate  functions  x(q,0)  and  y(q,0)  are  of 
the  same  order  of  magnitude  as  those  of  the  series  for  the  stream  function 
^(qj0)  (equation  (4l)).  By  analogy  with  the  preceding  section  the  first- 
order  terms  in  both  x  and  y  can  be  similarly  summed.  First,  let 
them  be  written,  for  q  <  1,  as 


2a*W”>(T)  cos  (n  —  1)0  —  1^^(T) 
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w 


cos  ne  sin  e  -  ein  n9  g,oa-£  j  (129) 


If  the  hypergeometric  functions  are  substituted  by  the  dominant  tern  of  their  asymptotic 
expansions  and  the  coefficients  are  approximated  by  then  x(q,6)  and  y(q. ,6)  reduce  to 
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Ao  sin  6  +  X(<it,e)  cos  9 


In  order  to  eliminate  the  error  introduced,  equations  (124)  and  (125)  can  similarly  be 
written  as 
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•where  Gn(t)  is  defined  hy  equation  (113)*  and  Oft, l(T )  is 
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the  following  equations  are  obtained  by  use  of  the  corresponding  func tions  xo(q„»0) 
y0(<l>e)>  and-  4>o(^0)  of  the  slmi:Lar  incompressible  flow: 


y(<i,e)  =  yi(<ij0)  +  ya(t^)  +  y3(<i>0)  (^3) 
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On  the  other  hand,  for  1  <  q.  <  V  hut  ------  <t  <  tjV  ,  the  dominant  terns  of  the 

£-P  +  JL 

bypergeometric  functions*  by  expressions  (226)*  (227)*  (233)*  (234)* 
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f(q,0)  =  [x(X,D  +  X(\,n)  +  2A0(jt  -  0)  J  cos  (ji°  +  +  |cr(X,l) 

—  a(X,Tj)l  sin  (\i°  +  ,  (¥?2) 
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Again  X2(q,0)  and.  y2(o.>®)  are  defined  by  equations  (144)  and  (145), 
but  the  differences  of  the  hypergecmetric  functions  involved  in  H±v(t) 
and  H^v^t)  are  now 

^±v(T)  =  F±v(t)  -  ^  T±V  cos  (voj  + 


^±v,l(T)  =  F±v,l(T)  -  4^  003  (v“  +  f) 

Improvement  of  Convergence  of  A^,  1^,  and 

As  seen  from  part  II,  as  soon  as  the  function  +o(q,0)  i3  specified, 

the  coefficients  of  the  inside  and  outside  series  can  be  chosen  so  as  to 
make  the  conditions  of  continuity  sufficient  for  the  sets  of  unknown 
constants,  for  example,  An  and  1^,  to  be  deteimined.  Inasmuch  as  the 

conditions  are  applied  at  the  circle  of  convergence,  it  is  generally  very 
tedious  to  evaluate  An  and  Dn  from  the  slowly  convergent  series 
(equations  (54)  and  (55))  when  Bn  and  Cn  are  given.  In  order  to  bring 
these  expressions  to  manageable  foims,  the  following  procedure  can  be 
used. 


2(3  +  1 


<  t  <  1 


xin 


By  considering  the  function  +o(q,0)  the  following  identities 
(appendix  F)  for  n  >  1  can  be  deduced  without  difficulty 


_  "*  XI  + 
n  0 

(153) 

1  i  00 

nAn  =  z  XI  — 

(154) 

0 

7  •  ;  t  p*  -  -  5 

(155) 

Furthermore,  the  sum  on  the  rightr-hand  side  of  the  second  expression  of 
equation  (52)  can  be  written  as 


WACA  TN  Wo.  1445 


51 


+3 

H 

a 

a 

t- 

J  *1/1 

p 

>'~c 

<  UJl 

J3 

1 

w 

H  O 

o 

ta  0 

0  ,£5 

U  — 
0  lf\ 


8LP 


«UJ° 

H|  t: 


H 

H 

O 

{u/t 

P 

n 

•H 

H 

ITS 

h 

H 

«w> 

vy 

1 

*d 

8UJ° 

H  |  ^ 


d. 

8U]° 

I — f)  ^ 

il 


in 

m 

ra 

ri 

g 

? 

*H 

P 

n 

§ 

•tf 

& 

«IlA 

ff) 

<1 

•> 

^J- 

in 

OJ 

H 

" — ■* 

n— ' 

P 

§ 

o 

*H 

P 

oj 

ro 

P 

& 

CVJ 

0 

w  d 

§ 

0 

U 

o 

•H 

0 

0 

• 

0 

OJ 

U  ' 

*N»vlt 

Pi 

H 

M 

'a 

0 

Ch 

0 

p 

in 

P 

U 

© 

0 

O 

£ 

nd 

P 

P 

o 

P 

EH 

t* 

O 

£n(Tl)  —  £-n(Tl)  An  —  —  2Do£-n(Tl)  “  an - 5 -  r  (^59) 


52 


MCA  TN  No.  1^5 


to  H 

•H  ©  -P 

sr  a 

O  'd 

■O  a  0) 

I  2  3 

•H  aJ 
©  CO  O 

A  © 

-P  O  ^ 
(0 

to  H  © 

♦H  cri  rQ 


w 

c 

*H 

•P 

•H 

& 

r'i 


2 

M 

+3 

u 

£ 


where 


NACA  TN  No.  1445 


53 


H  H 
H  H 

W  vT 


h  h 
t30  <H 


rH  ^ — s 

t-  H 

'h  i 

*  r> 
_>  PR 
P'H 


J 


<D 

60  03 


H 

O  V 


54 


MCA  TN  No.  1U45 


IV  -  APPLICATION  TO  CASE  OF  ELLIPTIC  CYLINDER 
The  Functions  zQ(w),  WQ(w),  and  A(w) 


An  irrotational  flow  of  an  incompressible  fluid  about  an  elliptic 
cylinder  with  a  circulation  ro  is  represented  by  the  complex  potential 

Wq(zo) • 


w0(z0)  =  5  + 


±  +  lr 


o 

2rt 


loge  £  with  zD  =  £  + 


d 


(1 64) 


where  the  flow  at  infinity  is  assumed  to  be  parallel  to  the  major  axis 
of  the  cylinder.  Here  all  the  quantities  have  been  rendered  dimensionless 
by  normalizing  t  by  a  length  a  and  W0(z0) ,  by  Ua.  Then  the  major 
and  minor  axes  are,  respectively,  1+e  and  1  -  e  ,  where  €  <1. 


By  differentiating  equation  (l64)  with  respect  to  zD,  the  dimension¬ 


less  complex  velocity  of  the  flow  is  obtained,  namely, 

i  - 1 


2  irQ 

5  +  0 


w  = 


2  7T 
£2-«2 


the  inverse  solution  of  which  is 


£(w) 


£o+  4(1-W)(1-  d*) 

2tt 


2(1  -  w) 


It 


(165) 


This  function  is  two-valued  with  two  simple  branch  points  at  w  =  1  +  o(tq2 J 
and  w  =  g  2  +  o(rQ2),  namely,  V  =  e-2  +  o(r q2).  (See  section  entitled 
"Analytic  Continuation  of  Solution  (Branch  Point  of  Order  1).")  The  principal 
value  may  be  defined  by  the  convention  that  — jt  <  arg  (1  —  w)  <  it 
and  1<  |w|  <  e”2  with  a  cut  Joining  the  -two  branch  pointB.  With  the 
principal  value  so  chosen,  the  domain  wherein  the  real  part  of  zQ  is 
less  than  or  equal  to  0,  excluding  the  interior  of  the  body,  corresponds 
uniquely  to  the  hodograph  given  by  W(zQ)  or,  conversely,  on  account  of 
symmetiy,  the  domain  wherein  the  real  part  of  zQ  is  greater  than  or 
equal  to  0  will  correspond  to  the  second  branch  of  the  function  £(w). 


By  substituting  £(w)  in  equations  (165)  and  (l6H),  an  expansion 
regarding  as  a  small  parameter  gives  in  accordance  with  equations  (35) 

4  Jt  n 

and  (36),  provided  that  ^  «  1, 
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Zo^(v)  =  - 


W0^°^(w)  = 


±jzs£)  1/2  +  c2  ( '±=±) 1,2 

1  -  w  /  VI  -  €  "/ 


^r-tetr 


zO(l)(w)  -  —i.'  1  e 


1  -  V  l  _  e2w. 


(166) 


•(167) 


(168) 


W0(l)(w)  =  -i 


1  -  w 
-  2ni 


—  +  logs  (1  “  w)  “ 


1  - 


-  logg  (1  -  e2v) 


(169) 


The  functions  z0(°) (w)  and  W0^°VW)  are  the  transition  function  and 
the  complex  potential,  respectively,  for  zero  circulation  (compare 

reference  1);  and  zQ(^V'w)  and  W0(^(w)  represent  the  firsi>-order 
contributions,  due  to  circulation,  to  z0(w)  and  W0(w),  respectively. 
On  separating  into  real  and  imaginary  parts,  it  was  found  that  from 

z0(°)(w)  and  WQ^°Vw)  there  result 


:o(°)(4,0)  »-i|[l(q,0)  +  J (<l,e)]l/2  +  e2[le  +  J  1(<l,©)]  1/2 


y0(o)(<i,e)  -  J 


and 


[-  1(4,0)  +  J(<i,e)]1//2  -  e2  [-I€  +  J  X(4,0)] 


1/2 


(170) 


(171) 


<p0(o)(q,e)  =  -  |[i(q,e)  +  J(i,e)]1/2  +  [i€  +  V2  j 

^0(o)(^e)  =  77=  {  [-  i(q»e)  +  J(q,e)]l/2  -  [-  i6(q,e)  +  J-1(<i,e)] 1/2 


(172) 


r  (173 ) 


where 


Kq,0) 


1 


(l  +  g  2)q,  cos  6  +  c242 
1  —  2q  cos  0  +  q_2 


(174) 
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Ke(q,0)  =  [(1  -  2^  oos  0  +  eS2)1/2  +  1  -  €2q  cos  ej  1 !2  (185) 

8!(<1,0)  =  Jjl  —  2q  cos  0  +  q.2)^2  —  1  +  q  cos  ©j^2  (186) 

K€(q,0)  =  [(1  -  2«2q  cos  0  +  e^q2)1^2  -  1  +  e2q  cos  0]  ^  (187) 


Expansions  of  W0(w)  and  zQ(w) 


As  the  domain  considered  is  complex  and  has  two  singularities  at 
v  s  1  and  v  =  the  function  WQ(  °)(w)  defined  by  equation  (167)  is 

single-valued  and  continuous  in  |w|  <  1  but  is  discontinuous  across  the 
cut  in  1  <  |v|  <  e~2.  The  expansions,  as  shown  in  reference  8,  are 

00 

W0(°)(v)  =  -  y~  A^w31  |w|<  1  (188) 

0 

and 

00 

W0*0)(v)  =  1  Y1  (Bnwv  +  Cnw“v)  1  <  |w  j  <  €~2  (189) 

0 

where  the  coefficients  are  all  real  and  are  defined  hy 

An  =  2SJ1)  -  (1  +  €2)Sn_1(i)  n  >  1  (190) 

Ao  =  2S0(1)  “  2 
Ajl  .=  0 

with 

s  w  ■  iyr(n-m  +  i)r^.ti)-  e2m  (191) 

n  *  0  r(n  -  m  +  l)r(m  +  1) 


2n 


and 


b,  =  KiW-(i*  *2)s„(°>]« 


092) 
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*o(o)(*,e)  »  (Bnq.V  +  Cnq.-^)  cos  ve  1<  q  <  6  2 


(200) 


Similarly,  from  W0(1)(w),  the  result  is 


00 

<PO^(<1,0)  =  -H  An^1^qn  Bin  nfl  -  2n 

2 

CO 

to^^(q,0)  =  —  y  Ap^^q11  cos  n0 

2 


q  <  1 


(201) 


(202) 


where  the  coefficient  An^^  is  given  "by 


2n 


(203) 


and. 


00 

<Po^1Hol,0)  =  ^  [Bn(l)qn  -  cn^<Tn]  s3-n  n6  -  («  +  0)  (204) 


o°  r 

+o(l'(a,e)  =51  [Bn(l)an  +  Cn(l)<l“n]  COS  n0  -  l0ge  q  +  1 


(205) 


where 


V 


(1)  .  (l  .  It31 


n 


(206) 


Cn(l)  =  1  +  ± 
1  n 


(207) 


The  expansion  of  zQ(w)  can  he  carried  out  in  like  manner,  but  it 
was  found  simpler  to  deduce  it  directly  from  equation  ( 37)  •  Then, 
corresponding  to  W0^°^(w)  defined  by  equations  (188)  and  (189), 
z0(°)(w)  is 


w  |  <  1  (208) 
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By  taking 


1 


2 


and  Mq_  =  0.60, 


Being  the  free— stream  Mach 


number,  the  compressible  flow  without  circulation  is  shown  in  figures  2 
and  3  in  both  the  t ,0—  and  the  x,y— planes.  The  flow  is  everywhere 
continuous  and  the  highest  Mach  number  attained  at  the  central  section 
is  1.24.  The  profile  calculated  in  this  case  is  symmetrical  but  is  far 
different  from  the  ellipse  from  which  it  was  derived.  The  thickness 
ratio  is  1:2  as  compared  with  3*5  of  the  original  elliptic  section.  The 
lateral  distortion  thus  is  much  more  pronounced  than  that  suffered  by  the 
longitudinal  dimension.  The  characteristic  feature  in  the  case  of 
transonic  flows  is  that  the  central  portion  of  the  derived  body  where  the 
flow  is  supersonic  is  always  flat  in  comparison  with  the  original  body. 


In  order  to  exhibit  the  ccmpressiblity  effect,  the  corresponding 
problem  of  the  incompressible  flour  over  the  same  body  must  be  solved. 
Technically,  this  does  not  offer  any  difficulty.  However,  in  order  to 
simplify  the  numerical  work,  an  incompressible  flow  resulting  from  a 
superposition  of  a  parallel  flow  on  a  source-sink  combination  has  been 
considered.  It  gives  rise  to  a  body  which,  by  adjusting  the  strength  of 
and  the  distance  between  the  source  and  sink,  approximates  closely  the 
given  body  with  an  error  of  about  2  or  3  percent.  (Compare  fig.  3-) 

The  pressure  distributions  over  the  same  body  for  both  compressible  and 
incompressible  fluids  are  compared  in  figure  4.  The  results  calculated 
according  to  the  Von  Kamnan-Tsien  and  the  Glauert-Prandtl  formulas  are 
also  given. 

The  wide  disagreement  between  the  exact  and  approximate  curves 
serves  as  a  proof  that  in  the  case  of  transonic  flows  the  thickness 
effect  of  the  body  is  no  longer  secondary  and  cannot  be  ignored  entirely, 
as  is  done  in  both  approximate  formulas.  The  immediate  cause  for  these 
discrepancies  seems  to  be  the  fact  that  the  point  of  equal  pressure  on 
the  surface  of  the  body  in  the  case  of  compressible  and  incompressible 
fluids  does  not  correspond  with  that  of  the  free— stream  pressure.  This 
(deviation  has,  in  fact,  been  observed  even  in  the  case  of  subsonic  flow 
about  a  Joukcwski  airfoil  of  thickness  ratio  O.lp  (reference  9)-  The 
reason  is  that  the  flow  of  a  compressible  fluid  over  a  body  is  partly 
compressive  and  partly  expansive.  By  increasing  the  free— stream  Mach 
number,  the  speed  of  the  flow  in  the  neighborhood  of  the  stagnation 
point  tends  to  decrease  and  the  speed  far  away  from  the  stagnation  point, 
that  is,  in  the  supersonic  region,  to  increase.  The  hodograph 
corresponding  to  the  zero  streamline  then  becomes  longer  and  flatter 
as  the  free— stream  Mach  number  increases.  Consequently,  the  point  of 
intersection  of  the  compressible  and  incompressible  holographs  will 
shift  toward  the  small  inclination  of  the  velocity  vector,  that  is, 
away  from  the  stagnation  point.  This  effect  will  be  more  pronounced 
for  larger  thickness  ratio  as  well  as  Mach  number. 

The  general  applicability  of  the  Von  Kaman-Tsien  formula  has  been 
questioned  and  examined  by  Tsien  and  A.  Fejer  (reference  10)  from 
purely  geometric  considerations,  and  the  limitations  therein  raised  seem 
to  be  further  substantiated  by  the  present  result.  The  question, 
however,  remains  to  be  answered  as  to  whether  critical  conditions  for 
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the  validity  of  the  Yon  Karman-Tsien  theory  could  he  given  in  terms  of 
the  thickness  ratio  and  Mach  number.  From  the  practical  point  of  view 
establishment  of  such  conditions  would  constitute  a  result  of  major 
importance. 

p 

For  the  case  of  circulatory  flow,  ^  is  taken  to  be  0.05.  Here, 

for  simplicity's  sake,  A^1)  =  1  and  B,-/^  =  (l  -  t^)-^  by  equation  (69), 

so  that  the  circulation  of  the  compressible  flow  is  (l  —  tj)  ^ro, 

according  to  equation  (74).  The  constant  C0^  could  have  been  determined 

by  successive  approximations,  but  in  this  particular  case  it  is  small 
enough  to  be  neglected.  The  calculated  flow  pattern  is  shown  in  figures  5 
and  6  for  both  planes.  The  body,  as  expected,  is  unsymmetrical  with  a 
negative  camber  of  about  2  percent.  The  highest  Mach  number  reached 
at  the  upper  surface  is  I.33  and  that  at  the  lower  surface  is  I.I5. 

The  lift  coefficient  is  O.65  corresponding  to  an  angle  of  attack 
at  zero  lift  of  about  2.2°.  For  camparsion  with  the  case  of  incompressible 
flow,  the  problem  is  simplified  by  considering  the  incompressible  flow 
over  the  body  resulting  from  the  superposition  of  a  small  negative  camber 
on  a  geometrically  similar  ellipse  of  thickness  ratio  1:2.  The  estimated 
lift  coefficient  of  such  a  system  is  approximately  0.45.  The  ratio  of 
the  two  coefficients  thus  is  nearly  1. 5,.  which  is  much  higher  than  the 
value  of  I.25  given  by  the  Glauert-Prandtl  formula  for  the  present  case. 
Actually,  if  the  oval  shape  is  considered,  the  incompressible  value  of 
the  lift  coefficient  would  have  been  lower  than  that  obtained,  owing  to 
the  fact  that  the  central  flat  portion  reduces  the  speed  of  the  flow  and 
hence  leads  to  a  higher  pressure.  Thus,  the  ratio  of  about  1.5  should 
be  regarded  as  a  lower  limit. 

Guggenheim  Aeronautical  Laboratory 

California  Institute  of  Technology 

Pasadina,  Calif.,  September  5*  19^6 
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APPENDIX  A 

PARTICULAR  INTEGRALS  OF  q>(q,e) 


Consider  the  total  differential  dq>  which  from  equations  (10 ) 
and  (11)  is 


dqp  =  -  -f  (1  -  M2)  dq  +  -f  d@ 


If  the  particular  integrals  of  if  are  given  hy  the  first  group  of 
expression  (24) ,  then 

a 


dq) 


_  _  v£o  f  _  y.2\  ♦v(a)]  -  sln  V0 
V  p  (1  M  )  q  1  COS  V0 

if 


dq 


P  dq 


nvF  (t)  I  i 
1  *V^J  jjsin 


cos  ve 

V0 


From  equation  (13)  ijr^(q)  can  he  eliminated  as  follows: 


dip  = 


Sq 


_o„V 


de 


Fv(t)|v(t) 


Vfv(t)iv(t) 


1 

r  ^ 

1 

!  sin  V6  l 

4 

—  cos  ve  f 

L  J 

dq 


sin  V0 
—  cos  V0 


de 


The  integration  therefore  gives 


A 


V  -  Vfv(t)|v(t) 


sin  VQ 
cos  V0 


>  +  Constant 


(220) 


Similarly,  the  following  equation  corresponding  to  the  second  group  of 
expression  (2i+)  is  obtained: 


Pn  -1/^  ,  v.  /  x  I  sin  ve 

qp  =  — q  (t  (t  )  {  —  cos  vQ 

P 


J 


>  +  Constant 


(221) 
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APPENDIX  B 
PROOF  OF  THEOREM  1 


Consider  a  small  neighborhood  of  q.  =  0  where  the  hypergeometric 
functions  can  be  accurately  represented  by 


Fn(T)  =  f(T)Tn(r) 

If  the  complex  potential  for  the  similar  incompressible  flew  is 


1  + 


■(1) 


ill 


n 


WD(w)  =  -  £ 
0 


w|  <  1 


let  W(w,t)  be  a  new  complex  potential  such  that  the  stream  function 
for  the  compressible  flow  is  given  by 

n<L,e)  =  3W  [W(W>T)] 

Then  W(w,  t)  must  be  of  the  foim 


W(v,T  )  =  _  ^>~  An(tw)n  +  [  f ^ 1  ^ (t  )  -  f^1^(Tl)]  —(^w)11 


>■  Iwl  <  1 


But  it  is  easily  seen  that 


in(t»)n  = 


Therefore,  W(w,t)  can  be  represented  by  two  closed  functions 


W(v,t)  = 


f(T) 

f(Tl) 


Wo(vt)  +  [f^(r)  -  f^^TijJjf  [w0(w)  +  A0]  ^ 


(222) 


The  function  represented  by  the  integral,  when  continued,  will  give  rise 
to  a  term  loge  tw.  No  longer  is  w  =  0  a  regular  point)  the  singularity 

must  be  either  at  w  =  1  or  at  origin. 
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APPENDIX  C 


ASYMPTOTIC  REPRESENTATION  OF  HYPERGEOMETRIC  FUNCTIONS 


The  asymptotic  integration  of  the  hypergeometric  equation  for  a  large 
positive  parameter  v  has  been  discussed  (reference  8).  It  was  shown 
that  its  solutions  under  this  condition  are  of  the  exponential  type  when  the 


variable  v  is  in  the  interval  8  —  r  —  - i — 

2p  +  1 


—  8,  where  8  >  0,  and  of 


the  oscillatory  type  when  t  is  in  the  interval  ~= — j  +  8  <  t  <1-  8. 

The  point  t  - 1 —  i3  a  singularity.  Thus,  to  the  first  approximation, 

2p  +  1 

the  hypergeometric  functions  for  the  interval  are 


Fv(t)  -  f(r)Tv(r)  J  1  +  o(i) 


F_v(r)  ~  f(r)T“  (t)  1 1  +  0(J 


v  >  N 


(223) 


(22b) 


where 


T(t  )  = 


f( t)  =  (l  —  r)a2A(l  - 

2  [ad-r^.d-aVt 

(1  +  tt)a(i  __  r)1/^  +  (1  —  a^r)1^ 


(225) 


Here  0^— J  in  each  case  denotes  the  fact  that  the  term  is  unifomly  of 
the  order  of  l/v  when  v  is  sufficiently  large  and  is  a  function  of  l/v 
For  the  interval  — - ~'-j-  <t  <  1,  the  hypergeometric  functions  are 


(226) 


F  (t)  f ( t )T  \ t )  cos  (V<0  -  y  )  1  +  01  — 


F_v(t)  ~  if(T)T  v(t)  cos  ^va>  +  1  +  of— 


(227) 


where 
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f(T)  =  2(1  -T)a2A(a2r  -  1)-1A 


T(t) 


2(gp)a^2  1 

(1  +  a)a  \j2pr 


> 


co(t) 


a  tan  ^  — 
a 


(228) 


(229) 


The  values  of  f(T),  T(t),  and  <d(t)  are  given  in  table  1.  Here 
the  parameter  v  may  he  any  positive  large  number.  When  it  is  a  large 
positive  integer,  these  formulas  will  automatically  represent  the 
hypergeametric  functions  defined  by  expression  (16). 

In  the  respective  domains  of  validity,  the  asymptotic  expansions 
may  be  differentiated  with  respect  to  t .  To  the  same  order  of  approxi¬ 
mation,  it  can  be  shown  that  for  V  >  N  in  the  interval 

Fv,i(T)  ~  s(T )tv(t )  [1  +  ofi)  j 

f-v,At)  ~  g(T),r_v(T)  j^i  + 


(230) 

(231) 


where 


6(t) 


2f  (t  ) 


(1  -  t)[1  +  |0(t)] 


and  in  the  interval 


2f3  +  1 


icA)  = 


<  t  <  1 


O 

1  -  a  t 


1  -  T 


Fv,l^T  )  ~  g(T  )tV^t  )  008  (V<J1  “  A-  f )  1  +  °(^) 


F_Vj1(t)  -  ig(r)T"V (t)  cos  (va>  +  |i°+  ^  1  +  °(y) 


(232) 


(233) 

(234) 
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where 


l(T)  = 


?2|3t  ( 1  —  t) 


(235) 


(!°(t)  =  COS 


-1.1  -  T 


Here  Fv  T )  and  F_V^^(T )  represent  respectively  the  functions 

F(av  +  1,  Tjv  +  1;  'cv  +  lj  t)  and  F(l  +  av  -  °v  +  I*  1  +  by  -  cv  +  lj 

2  —  cv  +  1;  t).  When  v  takes  integral  values  n,  expressions  (231) 
and  (234)  will  represent  asymptotically  F_n  ^(T  ) .  The  values  of  g(r) 

and  u°(t)  are  given  in  table  1. 

In  practice,  a  more  accurate  representation  for  the  functions  Fv(t) 
and  F_^(t)  and  their  derivatives  is  often  required,  especially  in  the 


interval  0  <  t  < 


The  fact  that  this  is  the  case  is  quite  evident 


2(5  +  1 

in  part  III.  In  order  to  cany  the  approximation  to  the  second  order, 
according  to  the  method  given  in  reference  8,  a  simple  evaluation  gives 


Fv(t)  ~  f(r )T  (t)  1  + 


•(l)(r) 


F  (t)  ~  f(T)T"v(T)  1  _  £ - ill  +  of— 

“V  V  Vv2 


0  <  T  < 


(236) 

. 1 . 

2p  +  1 

(237) 


where 


f(i)(T)  =  A,  [ <£i<L.  1  g)  (i_|0)  ioge  _(”ir- 

16  p  a  e  a  -  1 


i  (to-1  -  !)♦¥(!  -io-3) 


(238) 


Similarly,  by  differentiation  with  respect  to  t  of  expressions  (236) 
and  (237),  the  following  expressions  are  obtained: 
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Fv,l(T)  ~  8(T )Tv(t) 
r  ^t)  ~  s(t )t~v(t) 


(239) 


(240) 


* 


where 

g(D(T)  =  f(l)(T)  +  a!  (1  +  I0)(l  _  6c-2)  (241) 

By  comparison  of  the  values  of  f^d)  and  S^(t),  the  first  approxi¬ 
mation  of  Fv(t)  and  F_v(t)  is  seen  to  he  superior  to  that  of  Fv  j^(t) 

and  F  ,(t).  The  values  f(-^)(T)  and  6^)(t)  are  given  in  table  1. 

— Vj  X 

For  7  =  1.405  and  v  =  n  +  l/2  and  v  =  n,  n  being  a  positive 
integer,  the  two  groups  of  functions  Fv(t),  F_v(t)  and  Fv  ^(t), 

F_Vj1(t)  with  their  asymptotic  expressions  were  calculated  for  values 

of  t  varying  from  0  to  0.34  and  for  values  of  n  from  1  to  10.  The 
results  are  presented  in  table  2. 
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APPENDIX  D 

ASYMPTOTIC  EEPEESENTATION  FOE  |V(T )  AND  £_^,(t) 


Next  in  importance  are  the  functions  £v(r)  and  |_v(t),  defined  by 

equation  (29).  They  are  associated  with  the  particular  solutions  of 
cp(q,e),  and  eventually,  through  the  latter,  will  appear  in  the  various 
functions  in  the  problem  of  compressible  flow.  As  was  shown  in  part  II, 
the  whole  scheme  is  based  on  the  stream  function  i(q.,0),  and  its 
determination  depends  on  the  efficiency  of  the  determination  of  the 
coefficients  of  the  power  series  representing  'Kq*©).  In  order  to 
facilitate  such  evaluation,  the  asymptotic  expressions  are  again  powerful 
tools.  In  deriving  these  expansions,  it  is  also  convenient  to  start 
from  the  differential  equation  for  |v  (t  )  which,  from  equation  (13)., 
is  as  follows: 


C(0  + 


1  - 


tMt) 


JL 

2t 


CTO  -e 


]- 


(242) 


This  is  the  celebrated  Eiccati  equation.  This  proves  to  be  an  adequate 
form  for  the  asymptotic  development  of  !v(T)  in  "the  interval 

0  —  t  < - i — .  Suppose  the  expansions  are  of  the  following  forms: 

2P  +  1 


!v(t)  ~  lQ  + 


(243) 


I  (t) 


(244) 


where  |0(t)  is  defined  by  equation  (232).  Substituting  these  expressions 

in  equation  (242)  and  equating  the  coefficients  of  various  values  of  v-8 
to  zero  yields  for  the  interval 

|(1)(T)  =— (21,5) 

2(1  —  t)  '  / 

e(2)(t)  .  »2t(Bt  ±  2)  (1  _  t  -h\  _  S(1)(t)  (21,6) 

4|q(1-t)2V  '  2a2l0 

The  values  of  and  are  given  in  table  1. 
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APPENDIX  E 
PROOF  OF  THEOREM  2 


The  first  part  of  the  proof,  namely,  the  series 

*(<1,6)  =  21  sin  n0  q  <  1 

2 

vhich  is  absolutely  and  uniformly  convergent  in  any  closed  domain  in 
<1  <  1,  can  be  similarly  carried  out  (reference  1),  since  according 
to  equation  (159)  it  cun  be  deduced  that  for  large  values  of  n 

1^1 <  MM 

where  M  is  a  constant  independent  of  n. 1  Therefore,  it  follows  that 
for  tj_  «  —  ^  -  and  n  >  N 


|VnFn(r)M  <  M|V*l)n| 


For  the  series  (equation  (51))  it  is  noted  that  for  t  ^  +1 

and  v  >  N 


|FvC^)(x)|  <  Mtv 
F_v(r)(r)|  Mt—V 

where  the  region  - i - 5  <  ^  i -  +  S,  where  5  >  0  is  assumed 

2(3  +  1  2(3  +  1 

to  be  excluded.  Here  all  the  M's  ^are  different  but  independent  of  n. 
Furthermore,  from  equation  (l6o),  Dn  is  bounded  by 


Consequently,  for  n  >  N,  the  following  inequalities  are  obtained: 
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|Bn<iVFv(r)(T)|  <  M|Bn(<it)V 
Cn4”VBLv(p)(T)|<  M|cn(<itrV 

Pnq“^_n^(T)|  <  M  ' 


This  completes  the  proof. 
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APPENDIX  F 


DEDUCTION  OF  IDENTITIES  (153),  (15*0#  ^  (155) 


Consider  the  stream  function  of  the  similar  incompressible  flow 


+ota0)  =  Mn  sin  110 

2 


q  <  1 


tD(q,e)  =  ^  (Bnqv  +  Cnq[*j  cos  V0 


1  <  q  <  V 


At  the  circle  of  convergence  the  conditions  of  continuity  give 


*)T  sin  n0  =  (®n  +  cn)  cos 

2  0 


oo  00 

>  nAjj  sin  n0  =  y  v(Bn  —  Cn)  cos  v0 
2  0 

Since  the  limits  exist  by  hypothesis,  by  multiplying  both  sides  by  sin  n0, 
teim-by-term  integration  then  yields 


An  -  ^  XI  (®m  + 


(247) 


nAn  =  i  21  -  CjI. 

*  0 


(248) 


where 


I  =  — - —  +  — - — 
^  n  +  u  n  —  p 


u  =  m  +  — 

2 


By  use  of  appendix  B,  from 
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NACA  TN  No.  1445 


Figure  1.-  Mapping  of  half  plane  D. 
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Figure  2.-  Zero  streamline  of  a  compressible  flow  in  hodograph  plane.  €  =  Mi  =  0.60; 
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Figure  3.-  Comparison  of  derived  profile  for  €  =  —  and  =  0.60  with  profile  given  by 

Cl 

source -sink  combination  for  M.  =0. 


NACA  TN  No.  1445 


89 


Figure  4.-  Pressure  distribution  along  major  axis  of  ellipse. 
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NACA  TN  No.  1445 


Figure  5.-  Zero  streamline  of  a  compressible  flow  in  hodograph  plane.  €  =  — ;  =  0.60;  —  -  0.05. 
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